Can Modern Nuclear Hamiltonians Tolerate a Bound Tetraneutron? 
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I show that it does not seem possible to change modern nuclear Hamiltonians to bind a tetra- 
neutron without destroying many other successful predictions of those Hamiltonians. This means 
that, should a recent experimental claim of a bound tetraneutron be confirmed, our understanding 
of nuclear forces will have to be significantly changed. 1 also point out some errors in previous 
theoretical studies of this problem. 
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An experimental claim of the existence of a bound 
tetraneutron cluster (^n) was made last year Since 
then a number of theoretical attempts to obtain such 
bound systems have been made, with the conclusion that 
nuclear potentials do not bind four neutrons |M l^i 
However these studies have been made with simplified 
Hamiltonians and only approximate solutions of the four- 
neutron problem. In this paper I use modern realistic 
nuclear Hamiltonians that provide a good description 
of nuclei up to ^ = 10 and accurate Green's function 
Monte Carlo (GFMC) calculations to improve this sit- 
uation. (A list of earlier experimental and theoretical 
studies, also with generally negative conclusions, may be 
found in HHIi.) 

A series of papers 0, 0, 0] j have presented the develop- 
ment of GFMC for calculations of light nuclei (so far up 
to A = 10) using realistic two-nucleon (AW) and three- 
nucleon (NNN) potentials. For a given Hamiltonian, the 
method obtains ground and low-lying excited state ener- 
gies with an accuracy of 1 — 2%. I use this method in the 
present study; tests similar to those reported in the above 
papers have verified that the energies reported here have 
similar accuracies, with two exceptions: 1) when the en- 
ergies are very close to 0, the error is probably a few 100 
keV; and 2) the '^H calculations contain a technical diffi- 
culty that might be introducing systematic errors of up 
to 1 MeV. The reader is refered to |8j| for a review of the 
nuclear GFMC method or to the previously cited papers 
for complete details of how the present calculations were 
made. 

By using the Argonne vis NN potential (AV18) 9] and 
including two- and three-pion exchange NNN potentials, 
a series of model Hamiltonians (the Illinois models) were 
constructed that reproduce energies for A = 3 — 10 nuclei 
with rms errors of 0.6 — 1.0 MeV The best model, 

the AV18-hIllinois-2 (AV18/IL2) model, is used in the 
present study. 

GFMC starts with a trial wave function, '^t, which 
determines the quantum numbers of the state being com- 
puted. For p-shell nuclei studied in the above references, 
the Jastrow part of contains four nucleons with an 
alpha-particle wave function and A — 4 nucleons in p-shell 



orbitals. This is multiplied by a product of non-central 
two- and three-particle correlation operators. I use V&t 
for ^n with the same structure except there are two neu- 
trons in a ^5*0 configuration and two in the p-shell. The 
total J'^ of the '*n ground state is assumed to be 0+. 
There are two possible symmetry states in the p-shell us- 
ing LS coupling: ^S[22] and ^P[211]; both are used in 
these calculations. I could find no ^'t that gave a neg- 
ative energy for "^n using the AV18/IL2 model. GFMC 
calculations, using propagation to very large imaginary 
time (r = 1.6 MeV~^), also produced positive energies 
that steadily decreased as the rms radius of the system 
increased. 
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FIG. 1: Energies of *n in external wells versus the well-depth 
parameter, Vq. 



In a second study, I added artificial external wells of 
Woods-Saxon shape to the AV18/IL2 Hamiltonian and 
used GFMC to find the resulting total energies of the 
four neutrons. Figurenshows results for wells with radii 
R — 3, 6, and 9 fm (all have diffuseness parameters of 
0.65 fm) and varying depth parameter, Vq. It seems clear 
that four neutrons become unbound (have positive en- 
ergy) significantly before the well depth is reduced to 
zero. Linear fits to the least bound energies for each 
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Woods Saxon radius parameter are also shown; these ex- 
trapolate to an energy of +2 MeV when the external well 
is removed. (These least-bound solutions have large rms 
radii. A transition from the indicated linear behavior to a 
steeper linear behavior is observed for deeper wells; this 
transition is associated with a change to much smaller 
rms radii solutions. The steeper fits of course extrapo- 
late to much larger positive energies.) This suggests that 
there might be a ^n resonance near 2 MeV, but since 
the GFMC calculation with no external well shows no 
indication of stabilizing at that energy, the resonance, if 
it exists at all, must be very broad. In any case, the 
AV18/IL2 model does not produce a bound "^n. 

The authors of suggest that only small modifica- 
tions of existing nuclear Hamiltonians may be necessary 
to bind four neutrons. To study this possibility, I made 
a number of modifications to the AV18/IL2 model. In 
each case the modification was adjusted to bind "^n with 
an energy of approximately —0.5 MeV; the consequences 
of this modification for other nuclei were then computed. 
Four of the modifications are reported here: long- and 
moderate-range changes of the NN potential in the 
partial wave; introduction of an additional NNN poten- 
tial that acts only in total isospin T = | triples; and in- 
troduction of a NNNN potential that acts only in T = 2 
quadruples. In all cases the complete AV18/IL2 Hamil- 
tonian was used with the additional term. 

The strong-interaction part of the AVIS NN potential 
consists of one-pion exchange with the generally accepted 
value of /^/47r = 0.075, moderate-range terms that are 
associated with two-pion exchange but which have phe- 
nomenologically adjusted strengths, and a short-ranged 
completely phenomenological part. The potential is writ- 
ten in terms of operators which can be used to produce 
the potential for any partial wave. By making correlated 
changes to the radial parts of the central, cr^ • ctj, • tj, 
and (Tj • (TjTi ■ Tj operators, one can change the potential 
in just the S = 0,T = 1 partial waves: ^So,^D2,^G4, etc. 
A corresponding change to the four operator terms 
can cancel the change in the partial wave, leaving a 
change to basically only the ^Sq partial wave. 

In the first such modification of the AV18, I changed 
just the two-pion range part of the ^5*0 partial wave, so as 
to leave the theoretically well established one-pion part 
unaffected. Increasing this two-pion strength by 4.9% 
results in a ^n energy of —0.87(3) MeV. (The statisti- 
cal errors in Monte Carlo computed numbers are shown 
in parentheses only when they exceed unity in the last 
quoted digit.) As is shown by the points labeled "mod- 
%-27r" in Fig. [3 this changes the % phase shifts by 12° 
over a large energy range and produces a bound dineu- 
tron (the energy is —0.88 MeV). These changes far exceed 
those allowed by modern phase shift analysis. A some- 
what smaller change that produces a bound "^n can be 
made by using the AVI' potential 0| in the ^5*0 partial 
wave (and AV18 in the other partial waves); this results 
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FIG. 2: ^So phase shifts from AV18 and modifications to it. 
The lines show the pp, pn, and nn phase shifts for the un- 
modified AV18 while the symbols show the modified results. 



in a ^n energy of —0.52 MeV and about a 50% smaller 
change in the ^Sq phase shifts (the points labeled "mod- 
%-AVl'" in the figure). However again ^n is bound, 
this time with an energy of —0.42 MeV. Note that the 
one-pion-range part of the potential is also changed in 
mod-^So-AVl'. These ^n and ^n bound states are quite 
diffuse; the rms radii are respectively 2.8 and 3.6 fm for 
the two ^n cases and 7.3 and 10.3 fm for the '*n. The 
"^n pair distributions have a peak containing about two 
pairs with a structure close to that of the ^n pair distri- 
bution and a long tail. Thus the "^n looks like two widely 
separated dincutrons. 
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FIG. 3: Energies of nuclei and neutron clusters computed with 
the AV18/IL2 Hamiltonian with modified AW potentials (%- 
2tv and ^So-AVl') and with no modification (AV18), compared 
with experimental values for known nuclei. 

As noted, these modifications of the ^Sq potential to 
produce minimally bound tetraneutrons also produce 
dineutrons with about the same binding energies; thus 
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they are physically unacceptable modifications. Figure|3| 
shows that they also introduce large changes to the bind- 
ing energies of other nuclei; for example "^H is ~50% over- 
bound and is stable or almost stable against breakup 
into '^H-l-n-l-n as opposed to being a resonance in that 
channel Also six and eight neutrons form bound 

systems, although three and five do not.. 

The authors of Refs. 0, 1^.0 concluded that the non- 
realistic Volkov potentials [l^ do not bind ^n. However 
these potentials do have bound dineutrons and the above 
results suggest that they thus might bind ^n. To study 
this I made calculations using the first four Volkov poten- 
tials in all partial waves and no NNN potential. These 
potentials indeed do bind ^n with energies of —0.91, 
-1.04, -0.47 and -0.71 MeV while the ^n energies are 
-0.56, -0.60, -0.35, and -0.42 MeV, respectively The 
rms radii of the '^n systems are all about 11.5 fm, which 
may explain why these bound states were not discov- 
ered in Refs. |^ The variational energies for ^n 
with modifications to the AV18/IL2 Hamiltonian are pos- 
itive; that is, only with GFMC improvement does the 
system become bound. However for the simpler Volkov 
potentials, the ^'t already give negative energies and the 
GFMC just improves these energies. 

It must be emphasized that these bound ^n results do 
not at all support an experimentally bound "^n. The more 
than 35-year old Volkov potentials are not realistic; they 
produce bound ^n, with the same binding energies as 
their deuterons; they have no tensor or L-S terms; and 
they can not reproduce modern phase shift analyzes in 
any partial wave. The one thing in their favor is that, 
by having a space-exchange component, they introduce 
some saturation in p-shell nuclear binding energies; how- 
ever with just one radial form they are even simpler than 
the space-exchange AVX' introduced in Ref _llj. 

The above results show that it is not possible to bind ''n 
by modifying the ^Sq potential without severely disrupt- 
ing other nuclear properties. The next NN possibility is 
the '^Pj channel. The net effect of these is a small repul- 
sion in neutron systems. Setting this term to zero had 
very little effect on ^n; one would have to introduce sig- 
nificant attraction to bind ^n and then again many other 
nuclear properties would be unrealistically changed. 

It has been suggested that modifications to the NNN 
or NNNN potentials, which are experimentally much less 
constrained than the NN potential, could be used to bind 
^n. Timofeyuk added a central NNNN potential to bind 
^n, but found that it resulted in ^He being bound by 
about 100 MeV [3, However, as she suggests, one 
should try less disruptive things. A NNN potential that 
acts only in T — ^ triples would have the same effect on 
^n as one with no isospin dependence, but no effect on 
•^H and ^He because the contain only T = ^ triples. A 
NNNN T = 2 potential would also not effect ^He and 
6Li. 

To study such possibilities, I added potentials of the 



forms 

V^MT = |) = V, Y. [Y{n,)Y{r,,)]PiT=l), 

cyclic 

V,,ki{T = 2) ^ V^Y. [Y{n,)Y(r^k)Y{rki)]P{T^2), 

cyclic 

Y{r) = [l-e-(™''-) ]2, 

to the AV18/IL2 Hamiltonian. Here m.,^ is the pion mass, 
the P are projectors onto the indicated isospin states, and 
V3 and V4 were chosen to produce ''n with ~ —0.5 MeV 
energy. These forms have the longest range that is possi- 
ble from strong interactions; the cutoff makes the radial 
forms peak at 1.55 fm. Using more confined radial forms 
only increases the problems reported below. 

It turns out that the couplings must be quite large 
to produce the minimally bound ^n: V3 = —440 and 
V4 = —4750 MeV, which result in '*n energies of -0.60(5) 
and -0.55(6) MeV. This can be understood as follows. If 
the AW potential is used to bind ^n, the pairs can sequen- 
tially come close enough to feel the attraction; this allows 
the four neutrons to be in a diffuse, large radius, distribu- 
tion. However a NNN potential requires three neutrons 
to simultaneously be relatively close and thus the density 
of the system must be much higher. Indeed the rms radii 
of the *n for the Vijk{T = |) case is only 1.88 fm, while 
that for Vijki(T — 2) is 1.61 fm. Such small radii result 
in kinetic energies that are an order of magnitude more 
than those for the "^n systems bound by modified ^6*0 po- 
tentials; for the Vijk{T= |) case the expectation value of 
the kinetic energy is ~87 MeV, while those of the AW and 
NNN potentials are -49 and -38 MeV, respectively. (As is 
discussed in Ref. ,_5J , GFMC directly computes only {H) ; 
other expectation values involve extrapolations. Here I 
have reported the extrapolated potential values and sub- 
tracted these from {H) to get the kinetic energy.) 

The very large coupling constants for the Vijk{T = |) 
and Vijki{T = 2) potentials mean that they have a large, 
even catastrophic, effect on any nuclear system in which 
they can act. This is shown in Fig. 2) for example 
Vijk(T=^) doubles the binding energy of ^Li and triples 
that of ^He, while Vijki{T = 2), which can have no effect 
on ^Li, quadruples the binding energy of ^He. As noted 
before both of these potentials have no effect on *He. 
Both potentials make ^H stable by more than 25 MeV 
against •^H-|-n-|-n. However the most dramatic result of 
these potentials is that every investigated pure neutron 
system with A > 4 is extremely bound and in fact is 
the most stable "nucleus" of that A. For Vijk{T = |) 
the energies are -62, -220, and -650 MeV respectively for 
A = 5,6,8, while for V^jki{T = 2) they are -358, -1370, 
and -6690 MeV. 

These enormous bindings indicate that matter will col- 
lapse with such potentials. This is to be expected for 
purely attractive many-nucleon potentials. One should 
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FIG. 4: Energies of nuclei and neutron clusters computed 
with modified AW7V and NNNN potentials. 



add a shorter-ranged stronger repulsion to obtain satu- 
ration. Such a repulsion might improve the results for 
A > 6 nuclei. I attempted to study this by using a repul- 
sive term with Yukawa radial forms of range 2m,r- How- 
ever, in order to get any appreciable effect on ^He, the 
repulsive coupling has to be made quite large; this then 
requires at least a doubling of the attraction to still bind 
^n; this results in potentials that are so strong that the 
GFMC starts to become unreliable. The apparent im- 
possibility of correcting the A — 6 results by such a term 
may also be seen from the rms radii of the ^n reported 
above; they are smaller than the experimental value for 
^Li and reasonable ^He radii. Thus a short-ranged repul- 
sion that still leaves the "^n bound will certainly result in 
A — 6 nuclei with too small rms radii. 

In all of these cases I have made isospin-conserving 
modifications to the AV18/IL2 Hamiltonian; thus there 
have been T = 1 additions to the AW potential, or a 
T = I addition to the NNN potential, or a T = 2 ad- 
dition to the NNNN potential. It might be objected 
that the modifications should have been made only for 
nn pairs or nnn triples or nnnn quadruples since the 
nuclear force is least well determined for such systems. 
Such changes would mean much larger charge-symmetry 
breaking and charge-independence breaking potentials 
than are presently accepted. But even so, the changes 
to the AW force, if limited to just nn pairs, would still 
bind two neutrons which would change the experimen- 
tal scattering length from — 18 fm to a positive value. 



Such a nn potential would still bind ^n and ^n. I estimate 
that it would still increase the binding of by 3 MeV 
while it would have no effect on '^He. Thus the Nolen- 
Schiffer energy for the A = 3 system would be some five 
times too large. Many of the devastating effects shown in 
Fig. ^ would similarly persist even if the potentials were 
limited to nnn triples or nnnn quadruples. 

In conclusion, should the results of Ref. be con- 
firmed, our current very successful understanding of nu- 
clear forces would have to be severely modified in ways 
that, at least to me, are not at all obvious. 
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